Spectral Curves and Localization 
in Random Non-Hermitian Tridiagonal Matrices 
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Eigenvalues and eigenvectors of non-Hermitian tridiagonal periodic random matrices are studied 
by means of the Hatano-Nelson deformation. The deformed spectrum is annular-shaped, with inner 
radius measured by the complex Thouless formula. The inner bounding circle and the annular halo 
are stuctures that correspond to the two-arc and wings observed by Hatano and Nelson in deformed 
Hermitian models, and are explained in terms of localization of eigenstates via a spectral duality 
and the Argument principle. 



I. INTRODUCTION 



Hermitian tridiagonal random matrices are studied in 
great detail, and many results are available on spectral 
properties such as density, statistics and localization of 
eigenvectors. They appear in several models of physics, 
as Dyson's random chains, Anderson's models for trans- 
port in disordered potentials, Ising spin models with ran- 
dom couplings, /3-ensembles of tridiagonal random matri- 
ces. Hatano and Nelson^ introduced a beautiful method 
to study the localization of eigenvectors, by forcing an 
asymmetry of upper and lower nondiagonal elements. 
Then the eigenvalues are driven from the real axis to 
curves in the complex plane, in patterns that measure 
the localization length of the corresponding eigenvectors. 

Tridiagonal random matrices that are non-Hermitian 
from the start are less studied. They model systems with 
asymmetric hopping amplitudes 2 -"^, describe properties 
of ID random walks^ or the evolution of population 
biology . Their spectrum is complex. In this work we 
study how the Hatano-Nelson deformation modifies it, 
the occurrence of spectral curves and the connection with 
the localization of eigenvectors. 

Let us then consider complex tridiagonal matrices with 



M = 



ai b 1 

C2 '"• 

b n 



bn-l 



(1) 



where all matrix elements are i.i.d. complex random vari- 
ables. Here we use the uniform distribution in the unitary 
disk of the complex plane. This implies that the eigen- 
value density of the ensemble is only a function of the 
modulus of the eigenvalue. The eigenvalues of a sample 
matrix of size n = 800 are shown in FigJT] (left). 

We next consider two deformations of the matrix M , 




FIG. 1: Eigenvalues in the complex plane of a single non Her- 
mitian tridiagonal matrix of size n = 800. The same matrix 
entries {at, bk,Ck} are used, with £ = and ip — (left), and 
£ = 0.5 and ip = (right). Note that the eigenvalues out of 
the circle have the same positions in the both cases. 



by a complex parameter z = e^ +l(p : 

ai 6i z n ci 



M{z n ) 



C-2 



M b (z) 



bn/z n 

ai b\/z 



ZC 2 



b n /z 



b n -l 

ZC\ 



bn-l/z 

a„ 



(2) 



(3) 



The two matrices are similar, Mj,(z) = SM(z n )S^ 1 , 
through a diagonal matrix with entries Sa — z l . The 
balanced matrix Mf,(z) is more convenient for numerical 
work. Since the matrices share the same set of eigenval- 
ues, a rotation of z by 2ir jn does not change the eigen- 
values of Mb(z). 

The eigenvalues of Mb(z) are shown in FigQ] (right). 
The distribution looks remarkable: a "circle" centered in 
the origin bounds an outer annular halo where the eigen- 
values appear exactly in the same positions as those in 
the left figure. The inner region is void: all the eigen- 
values that were there before deformation (£ = 0) have 
moved to the boundary circle. As \z\ becomes larger, see 
FigEl the circle enlarges as well, but the eigenvalues in 



FIG. 2: The same random matrix entries as in Fig[T] with 
£ = 0.7 and tp = (left), £ = 1 and tp = (right). 




FIG. 3: Motion in the complex plane of the eigenvalues of 
a single non Hermitian tridiagonal matrix of size n — 100 
as the phase tp is changed in half of its range. The same 
matrix entries are used, with £ = 0.5 (left) and f = 1 (right). 
Eigenvalues not belonging to the loop are seen to be fixed, 
corresponding to localized states. The loop becomes closed if 
the whole angular range is evaluated. 



the annular halo do not move, until they are swept by the 
circle. For large £ only the circle remains. This is not sur- 
prising: in the limit of large \z\ the matrix Mf,(z) simpli- 
fies to bidiagonal. The eigenvalue equation can be solved 
explicitly and gives E n — z n b\ ■•■&„. Then the eigenval- 
ues Ek = |z|e( ln ' 6 ') expi(0 + 2-Kk/n) are equally spaced 
and lie on a circle of radius r such that log r — £+ (log | b | } . 

Eigenvalues on the circle and in the halo respond differ- 
ently to the phase tp = argz. As tp sweeps the Brillouin 
zone from to 2ir/n, only the eigenvalues sitting on the 
circle move (and remain therein), while the outer ones do 
not have measurable changes at all. This is illustrated in 
Figj3l which also shows that an eigenvalue on the cir- 
cle moves to the position of a neighboring one as tp is 
increased by 2ir/n. By increasing the size n of the ma- 
trices, the "circle" is seen to become independent of the 
sample and more regular, (Fig[4]). 

The phenomenon described is analogous to what 
Hatano and NelsonA discovered for random tridiagonal 
Hermitian matrices (afe real, Ck+i = &|) where the unde- 
formed eigenvalues (£ = 0) are real. The deformation 
forces them to move into the complex plane and dis- 
tribute along a two-arc loop, with possible external wings 
of unaffected eigenvalues in the real axis (Figj5|) . The 
two-arc loop and wings of the Hermitian model corre- 
spond to the circle and annular halo of the non- Hermitian 
model discussed here. 



FIG. 4: Eigenvalues in the complex plane of 3 different sam- 
ples of size n = 400 (left) and n = 1200 (right). All matrices 
with £ = 0.5 and tp = 0. 




FIG. 5: Eigenvalues of a single Hatano-Nelson tridiagonal ma- 
trix of size n = 600, with £ = 1, bk = c& = 1 and random 
numbers at uniformly distributed in (—3.5,3.5). They be- 
long to two arcs or two wings in the real axis, which are the 
residuate of the undeformed spectrum. 



In Section II and III we study the spectral density of 
the undeformed ensemble and the localization of eigen- 
vectors, measured by the Lyapunov exponent or by the 
variance. In Section IV we explain the observed spec- 
tral features of the deformed ensemble by means of the 
Argument Principle of complex analysis and a spectral 
duality between the eigenvalues of M(z n ) and those of 
the transfer matrix. 
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FIG. 6: Eigenvalue density pod^l) of the undeformed ensem- 
ble (matrix size n = 1000, 100 samples). The dark line is 
obtained by smoothing the hystogram data on local windows 
of 10 bins (out of 300). The plateau is fitted by the value 
0.193(2). 



FIG. 7: The Lyapunov exponent 7(|-E|) evaluated numeri- 
cally (100 matrices of size n = 1000). The two curves are 
a quadratic fit near the origin, 0.2914(1) + 0.309(0) 2 , and 
the exact analytic expression lg(|£|) + 1/2 of 7 for \E\ > 2. 
The fit near the origin is consistent with the value found for 
po(0) (via Poisson equation). 



FIG. 8: Pairs (\E a \,y a ), where E a are eigenvalues and -y a 
are rates of exponential localization, for the eigenvectors of a 
matrix M of size n = 800 (£ = 0). The continuous line is the 
Lyapunov exponent 7(|_E|). 



II. THE SPECTRUM OF M 

Since the matrix entries of M are chosen to be uni- 
formly distributed in the unit complex disk, the av- 
erage eigenvalue density of the ensemble, po(E) = 
(- J2i feiE — Ei)), depends on \E\. The eigenvalue equa- 
tion 



Ck u k -i +a k u k + b k u k 4 



Eu k 



(4) 



written for the component Uk with highest absolute value, 
implies the inequality | E \ < |afc| + |6fe| + |cfc| < 3. Thus the 
disk that supports the density has a radius not exceeding 
3; the numerical evidence is that it has length 2. 

We diagonalized 100 matrices of size n = 1000 to ob- 
tain numerically the density of eigenvalues po(l-El) shown 
in Figj6] The lowest moments pk = (-kJ2i\Ei\ k ) = 
2-k Jq 2 dxx k+1 po(x) are also evaluated: p\ = 0.9107(12), 
A* 2 = 0.9678(22), p 3 = 1.1327(36) and p 4 = 1.4204(59). 



III. THE LYAPUNOV EXPONENT 

The numerical evaluation of the eigenvectors of M 
shows that they are strongly localized for all eigenvalues; 
indeed, for large matrix size, they decay exponentially 
(Anderson localization). The rate of decay is measured 
by the Lyapunov exponent, an asymptotic property of 
transfer matrices. 

The transfer matrix of a realization M is the product 
of 2 x 2 random matrices: 



HE) = n 



b^(E-a k ) -b^c k 



Its eigenvalues can be written as z\ 



large n the exponents £,±(E) become opposite: £ + 

= iiog|dett(#)| = iELii°sM-iogN^o. 



(5) 
For 



According to the theory of random matrix products, 
for large n the positive exponent £ + becomes independent 
of n and the realization of randomness, and converges 
to the Lyapunov exponent of the matrix ensemble. The 
Lyapunov exponent can be evaluated by an extension of 
Thouless formula to non-Hermitian matrices^, 



J <PE'p {E')\og\E-E'\ 



<log |6|>, (6) 



where po is the eigenvalue density of the ensemble of ma- 
trices M. Note that for complex spectra, the equation 
for 7 implies the Poisson equation V 2 7(i?) = 2npo(E). 
Therefore j(E) can be understood as the electrostatic 
potential generated by a charge distribution in the plane 
with density po(E). 

For a distribution of matrix entries that is uniform in 
the unit disk, it is (log|6|) = —1/2 and po is rotation 
invariant. Then the integral can be simplified: 

7 (|£|)= log|£|A/- (|£|) 

/'OC 
dE'E'po(E') log E' + 1/2. (7) 
,E\ 

The integral log |r — r'e lv \dip — 27rlogmax(r, r') was 
used. Afo(\E\) is the fraction of the spectrum inside the 
disk of radius \E\. For \E\ larger than the spectral radius 
it is 



7 (|£|) = log(|£?|) + 1/2, \E\>2. 



(8) 



The Lyapunov exponent is an increasing function of 
\E\. Its numerical evaluation is shown in Fig [71 

We checked numerically the exponential decay of eigen- 
vectors with a rate given hy j(\E\). If u is an eigenvec- 
tor of Af, with components {wfc}£ =1 , the numbers |ufe| 2 
provide the probability distribution for the position of a 
particle in the lattice 1 . . . n. We choose to measure the 



4 



400 r 
350 - 

300 - 




FIG. 9: Variance of eigenvectors (z-axis) and corresponding 
eigenvalues [x = ReE, y — ImE) of a random matrix Mj,(e^) 
of size n = 800 and £ = log 2. The radius of the hole is 
approximately 1.16. 



localization of the eigenvector by the variance in position: 
( n - V /2 

var ^ = \y^\u k \ 2 (k-k) 2 \ , (9) 

where k = J2k k\ u k\ 2 is the mean position of the particle. 
Other measures could be used, as the inverse participa- 
tion ratio or the Shannon entropy; in this case the eigen- 
vectors are peaked on small intervals, and the variance 
has a clear meaning. 

For an ideal state v that is exponentially localized, 
\vk\ 2 = (tanh.7) e 2 " y 1^1, (777, 3> 1), the variance is 
varju] = sinh(l/7). We use the same relation to com- 
pute a rate j a from the numerically evaluated variance 
of an eigenstate u a . In FigJS] we plot the numerical pairs 
(|E a |,7 a ) for the eigenvalues and eigenvectors of a single 
matrix M of size n = 800, together with the Lyapunov 
exponent 7(|.E I |), given by Thouless formula |J7J|. The 
numerical data are consistent with the picture of expo- 
nential localization of eigenvectors. 



responding complex eigenvalues (horizontal plane). The 
boundary of the circular hole is populated by the eigen- 
vectors which are delocalized. 

The existence of an empty disk and a halo of fixed 
eigenvalues for the deformed ensemble thus reflects the 
localization properties of the eigenvectors of M as a func- 
tion of \E\, i.e. the function 7(|_E|). 

Proposition: in the large n limit, Mb(e^ +lLp ) has no 
eigenvalues in the disk of radius r, where 

7 (r) = i (10) 

Proof: The hole in the spectrum of Mj,(z) can be un- 
derstood via the Argument Principle of complex analy- 
sis: the number of zeros of the analytic function f(E) = 
det[E — Mb(z)] inside a disk of radius r is equal to the 
variation of &rgf(E) /2ir along the contour of the disk. 

The function f(E) is related to the eigenvalues z±(E) 
of the transfer matrix t(E) by a duality identity^^: 

det[E - M b (z)} = •••&«) det[t(£) - z n ] (11) 

z n 

Then 

argdet[£ - M b (e i+tlp )} = const. 
+arg[e" ( «+-« )+m(v +- ¥,) - 1] 
+axg[e n ^—^ +i ^—^ - 1] 

Let us fix £ > and take the large n limit. Then: £ + > 
and < 0; &r g[e n ^+^^ +m ^+-^ - 1] equals nUp + - <p) 
if £+ > £, and n if £ > £ + ; avg [ e ^~-i)+^{v-v) - 1] = 
7r always. We also identify £+(E) with the Lyapunov 
exponent 7(|£'|). The variation of a,rgf(E)/2ir along a 
circumference of radius r is zero if £ > j(r). □ 

Since 7(0) is nonzero, there is a threshold value £ m .;„ ~ 
0.291 below which no hole opens in the spectral support. 

V. CONCLUSIONS 



IV. HOLE, HALO AND LOCALIZATION 

When the parameter £ is switched on, it appears in the 
corners of the matrix M(z n ) and modifies the boundary 
conditions in Q. In the transition from M to M(z n ), one 
expects that the eigenvalues of enough localized eigen- 
states do not change. An eigenstate u — {u^} of M(z n ) 
corresponds to an eigenstate Su of Mb(z), with compo- 
nents z k Uk- If \v,k\ ~ e for large k, the factor e k ^ 
delocalizes it if 7 < £. This simple argument by Hatano 
and Nelson indicates a threshold value 7(|-E|) = £ at 
which eigenvalues must be drastically influenced by the 
deformation. 

In Fig[9]we plot the variances (z axis) of the eigenvec- 
tors of a matrix Mf,(z — 2) of size n = 800, and the cor- 



The Hatano Nelson deformation opens a hole in the 
spectrum of non Hermitian tridiagonal random matrices 
with i.i.d. matrix elements. The eigenvalues that are 
swept to the boundary of the hole correspond to states 
that are no longer Anderson localized. This is explained 
in terms of a spectral duality, stability of the Lyapunov 
exponent, and the Argument Principle. 

Tridiagonal matrices with different strengths of ran- 
domness in the three diagonals would also show similar 
spectral features. 
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